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Cavity-assisted storage and retrieval of single-photon wave packets in optically thin spatially 
extended resonant materials are analyzed. It is shown that the use of cavity tuning allows one 
to store and recall time-symmetric double-sided exponential pulses with near unit efficiency. The 
optimal regime of the cavity tuning is determined and the effect of time jitter on the storage efficiency 
is analyzed. 
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I. INTRODUCTION 

A quantum memory for light may be defined as a con- 
trollable delay line for quantum states of the electromag- 
netic field. Such devices are currently considered as a ba- 
sic ingredient for scalable all-optical quantum computers 
[1] and quantum repeaters [2] , which makes them a topic 
of active research. Although the optical quantum stor- 
age can be realized by means of usual delay lines [3, 4], 
the most promising approaches to the problem involve 
interaction of light with an atomic system controlled by 
additional fields or pulses [5-13]. An important class of 
optical quantum memories is designed for single-photon 
quantum states which can be stored and recall using in- 
teraction with a single atom [6] and with an ensemble of 
atoms in the regime of electromagnetically-induced trans- 
parency [7, 8, 14, 15], stimulated Raman absorption [13] 
and photon echo [5, 10, 16-22]. Regardless of the ap- 
proach, the possibility of manipulating the atomic sys- 
tem strongly depends on its physical properties, which, 
generally speaking, restricts the set of materials suited 
for quantum storage. For example, the photon echo gen- 
erated using controlled reversible inhomogeneous broad- 
ening in a solid-state material [16-20] is of particular in- 
terest since it allows storage and retrieval without ap- 
plication of additional laser pulses. But the frequency 
manipulation via the linear Stark effect is possible only 
if an optical center occupies a non-centrosymmetric site 
in a medium (see the review in [23]). Besides, a signif- 
icant difference between the excited- and ground-state 
dipole moments is desirable. In this respect the scheme 
for quantum storage on subradiant states in a solid-state 
material [12] may be of interest. The advantage of the 
scheme is that storage and retrieval of a single-photon 
state may in principle be implemented without any ad- 
ditional control fields or pulses acting on the resonant 
atomic system, which makes it possible to use materials 
that are difficult to control. In the present work we fol- 
low this approach and consider the possibility of such a 
passive manipulation of atomic states in a tunable cavity. 
Enclosing an atomic ensemble in a cavity makes it possi- 
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ble to achieve high efficiency of quantum storage with op- 
tically thin materials, thereby reducing phase relaxation 
rates and parasitic absorption, to the extend possible. 

Atom-cavity systems have been investigated in the con- 
text of quantum memory in a series of papers [6, 7, 24- 
29]. In particular, it was shown that time-symmetric 
single-photon wave packets can be stored and recalled 
with nearly unit probability, provided that atom-cavity 
interaction is controlled by special laser pulses [6]. Fur- 
thermore, the storage and retrieval of photons in a A- 
type atomic medium enclosed in a running-wave cavity 
were optimized for different approaches to quantum stor- 
age [28]. It was shown that the optimal storage process 
is the time-reverse of retrieval and the highest possible 
total efficiency of a quantum memory device is equal 
to (C/(l + C)) 2 , where C is the cooperativity param- 
eter. The same efficiency can be achieved for the squeez- 
ing transfer between an atomic ensemble and fields in 
an optical cavity [24]. Besides, the storage efficiency 
was analytically studied for two-level and various three- 
level atomic configurations and for various single-photon 
pulse shapes [27]. In the present work the capabilities 
of the storage and retrieval via cavity tuning are consid- 
ered. The simplest strategy for light storage and retrieval 
involves steplike modulation of the atom-field coupling 
strength via cavity tuning [29]: the coupling is switched 
off after absorption of a light pulse and switched on be- 
fore the emission. In this case, the optimal pulse shape 
leading to the highest possible efficiency is asymmetric 
one and the output pulse is the time-reversed replica of 
the input pulse [28, 29]. In this paper, a simple scheme 
for manipulation of cavity parameters is developed such 
that symmetric double-sided exponential pulses can be 
stored and recall with a near maximum efficiency. The 
cavity tuning is apparently the counterpart of Stark shift- 
ing in an atomic system. In this connection it should be 
noted that in Ref. [30] Stark shifting in an atomlike sys- 
tem was used to control its interaction with emf in a 
cavity and possibilities of pulse shaping during sponta- 
neous emission were investigated by quantum trajectory 
simulations. In comparison with that work, we consider a 
full storage cycle — i.e., the storage followed by retrieval 
— and the optimal cavity control during the storage and 
retrieval of information is described analytically without 
the use of adaptive algorithms. 
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The paper is organized as follows. In Sec. II, we intro- 
duce a physical model of cavity-assisted interaction be- 
tween an extended atomic ensemble and a single-photon 
wave packet and present a general solution of the evo- 
lution equations in the case of dynamic cavity tuning. 
In Sec. Ill, the optimal regime of the cavity tuning that 
allows one to store and recall double-sided exponential 
light pulses is described analytically, the efficiency of the 
storage followed by retrieval is calculated and the effect 
of time jitter on the storage efficiency is analyzed. 



II. THE MODEL AND BASIC EQUATIONS 

Consider a system of N ^> 1 identical two-level atoms 
which are placed in a single-ended ring cavity and inter- 
act among themselves and with the external world only 
through the electromagnetic field (Fig. 1). We are inter- 
ested in the interaction of the atomic system with a reso- 
nant single-photon wave packet propagating through the 
cavity. In doing so we suppose that the bandwidth of col- 
lective spontaneous emission in the cavity as well as that 
of the incoming light are much smaller than the cavity 
lincwidth. Such a cavity-assisted interaction corresponds 
to the strong "bad cavity" limit [31] and can be con- 
sidered in the Born-Markov approximation. Besides, we 
restrict the consideration to a single cavity mode which 
subtends a small solid angle at the atoms, a cav , so that 
spontaneous emission rate into free-space modes is actu- 
ally not modified by the cavity. On the other hand, a cav 
is supposed to be larger than the diffraction angle of col- 
lective atomic emission so that the cavity may be called 
"bad" in the angular domain as well. 

We leave out any internal losses inside the cavity and 
consider only those through the partially transmitting 
mirror. The corresponding input-output relations for the 
field creation operators, provided that the mode wave- 
vector k is in a cav , may be written as 
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where g(kp) — rcxp(-ikp), r = yl — t 2 , t is the ampli- 
tude transmission coefficient of the semitransparent mir- 
ror and p is the optical path length for one round trip 
through the cavity, which takes into account all the re- 
fractive indices inside the cavity. The operator 5,J av is 
not normalized as a creation operator, but includes the 
enhancement of the cavity field and its phase shift with 
respect to the input field. As a result, the relations (1) 
and (2) have the same form as those for the complex 
field amplitude (see, e.g., Refs. [32, 33] for details). Here 
a| av is defined for the central point of the atomic ensem- 
ble which is separated from the coupling mirror by the 
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where T = W\/r/(l — r) is the cavity finesse and it is 
assumed that T ~> 1 . 

In the "bad cavity" limit, the values of the functions 
C{kp) and fi{kp) may be considered as constants within 
the atomic resonant line, which vary in accordance with 
Eqs. (3) and (4) under changes of p in time. It is con- 
venient to introduce the notations C(t) = £(kp(tj) and 
f3(i) = fi{kp{t)). The value of £(t) determines the atom- 
field coupling strength, while that of (3(t) describes phase 
shifts between the input and output fields. Therefore, in 
order to simplify the interaction and avoid unnecessary 
phase modulation of input and output photons, we put 
phase modulators PM3 and PM4 at the input and out- 
put of the cavity, respectively, each create phase shift 
exp(— if}{t)/2) upon the symmetrical modulation of the 
cavity optical length by modulators PM1 and PM2. As 
a result, the input-output relations become 
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Let us denote the ground and excited states of jth 
atom (j = 1,...,N) by |0j) and |lj), its position by 
rj and resonance frequency by uo- We assume that (i) 
the atoms are not moving as, for example, impurities 
embedded in a solid state material, (ii) the excitation 
volume may be approximated by a cylinder with the cross 
section S and the length L, and the Fresnel number of 
the excitation volume F = S(LX)^ 1 > 1, (iii) the atomic 
system has a small optical thickness, aL <C 1, where a 
is a resonant absorption coefficient. In the "bad cavity" 
limit, the Hamiltonian of the system, in the interaction 
picture and rotating-wave approximation, is 

H = ^..V^W&Jok,, ^ e *(«o-<")* + H.c. (7) 
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Here <?k,s = \ (^t 7 ) ' £k >s) is the atom-field cou- 

pling constant, bj = \0j)(lj\ is the atomic transition op- 
erator, flk, s is the photon annihilation operator in the 
radiation field mode with the frequency w = kc and po- 
larization unit vector e^ s (s = 1,2), V is the quantiza- 
tion volume of the radiation field (we take V much larger 
than the volume of the atom-cavity system), d is the 
dipole moment of the atomic transition, and 



£k(t) 



C(t), k in a ce 
1, otherwise. 



(8) 



For simplicity we assume that the vectors Sk, s and d are 
real. 
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FIG. 1: (Color online) Schematic of a quantum memory de- 
vice with a tunable cavity. Ml and M2 are perfectly reflect- 
ing mirrors, M3 is a partially transmitting mirror. The op- 
tical path length of the cavity is controlled symmetrically by 
phase modulators PM1 and PM2, while modulators PM3 and 
PM4 are used for controlling input-output phase relations. All 
phase shifts may be implemented by movable mirrors thereby 
reducing the losses. 



We define the following collective atomic operators: 

N 



-ikn-r,- 
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where ko is directed along the cavity mode within the 
sample, |ko| = Wo/ c j so that the general form of the state 
of the system can be written as 

|V(t)) = E/k,aW|0)|lk )S ) + P(t)|l}|vac) (10) 



with normalization condition ^ k \fu,s(t)\ 2 + \P(t)\ 2 = 
1, where |0) = |0i, 0a, . . . , On) is the ground state of the 
atomic system, |vac) is the vacuum state of the radiation 
field, |1 M ) = 4 iS |vac) and |1) = N-^R^O). 

Substituting Eqs. (7) and (10) in the Schrodinger equa- 
tion we obtain 



dt 

8P(t) 
dt 
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where </>(x) = N~ x ^\ exp(ix ■ r 3 ) is the diffraction func- 
tion. 

The photon density for the incoming wave packet at 
the input of the system reads 



(13) 



and for the emitted radiation we have the analogous 
equation with F ln (t) and /k,s( — °°) replaced by F out (t) 
and fk,s(t), respectively. Then the solution of Eqs. (11) 
and (12) may be written as (see Appendix) 



F out {t) = F in (t) + y/2-yC(t)P(t), (14) 



P{t) = P(-oo) exp (- J 7(C(r)+l)dr 

^C{r)F- m {r) exp (- j 7 (C(r') + 1) dr') , 



where C(t) = C(t)Nfi and 7 = 1/T 2 = l/2Ti. 
/i = 3\ 2 (8irS)~ 1 is a geometrical factor [34], and 
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Here 



(16) 



is the Einstein A coefficient. The function C(t), which 
can be expressed as C(t) — C(t)aL/A, is so-called co- 
operativity parameter and represents the effective opti- 
cal depth of the medium in the cavity. The maximum 
value of the cooperativity parameter, C m — aLT /2ir, 
is attained at exact resonance and is equal to the op- 
tical length of the medium multiplied by the number 
of round trips during the photon life-time in the cav- 
ity. The strong "bad cavity" limit is now defined as 
Cm7 *C Aw cav = 2irc/J : 'p. According to the definition 
of C(t) and Eqs. (3) and (4), there is a one-to-one corre- 
spondence between the values of C(t) and those of p(i) 
and (3{t) on each side of the cavity line. Thus, the cavity 
tuning is fully specified by the function C(t). 



III. STORAGE AND RETRIEVAL OF 
DOUBLE-SIDED EXPONENTIAL PULSES 

In this section we show how to tune the cavity length 
for the storage and retrieval of symmetric double-sided 
exponential pulses and calculate the efficiency of such 
a cavity-assisted quantum memory. There are several 
reasons for choosing such pulses for the analysis: first, 
symmetric pulses are suited for multiple storage and re- 
trieval without resorting to additional time reversal; sec- 
ond, single-photon pulses of such a shape may be gener- 
ated via spontaneous parametric down conversion when a 
nonlinear crystal is enclosed in a cavity and continuously 
pumped [35, 36]; third, such pulses are near optimal for 
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quantum-information processing from the viewpoint of 
tolerance against the effect of time jitter [37]. 

To begin with, we determine the time dependence of 
the cooperative parameter C(t) which leads to the spon- 
taneous emission of the double-sided exponential pulses 
(Fig. 2). Suppose that the spontaneous emission begins 
at the moment of time t = when atoms are coherently 
excited so that -P(Q) = 1. Let the maximum of the pulse 
amplitude falls at time t = T when the cooperativity 
parameter reaches its maximum value C m by the cavity 
tuning and the latter does not change after that. In this 
case the output pulse can be written as 

F out {t) = F m exp{-j{C m + l)\t-T\), t>0. (17) 

The trailing part of the pulse (t > T) is the usual ex- 
ponential decay when the cavity is tuned to resonance 
with the atomic transition. The leading part is the time- 
reversed exponential decay which the dependence C(t) 
should be determined for. For this purpose, according to 
Eqs. (14) and (15), we need to solve the equation 



V2 7 C(i) exp 



7(C(r) + 1) dr 



F m exp(7(C m + l)(t-r)), (18) 



provided that < t < T. By differentiating the loga- 
rithm of Eq. (18) we obtain Bernoulli equation for C(i), 
which has the following solution: 



C(t) = 



C n 



A{C m + 2) e -27(c m +2)t _ 1 ■ 



(19) 



The constant A is found from the condition C(T) = C„ 
to be 



A = e 2 7 (C m +2)T ( 1 + J_ 
Cm + 2 C,r. 



(20) 



Finally, substituting Eqs. (20) and (19) into the left side 
of Eq.' (18) we find 
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(21) 

Thus, if T = 0, then F m — \J2^C m and we have an 
ordinary exponential decay pulse with the highest pos- 
sible amplitude. On the other hand, when ~fC m T 3> 1 
then F m = ^/-iC m (C m + 2)/(C m + 1) exp(-T/T 2 ) and 
we have a double-sided exponential pulse the maximum 
amplitude of which decreases with T due to the phase re- 
laxation caused by spontaneous emission into free-space 
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FIG. 2: (Color online) The input, F in (t), and output, F out {t), 
double-sided exponential pulses (above) and corresponding 
cooperative parameter C(t) (below) as functions of time. 



modes. The norm of the output pulse is 
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and to the lowest order on 1/C m we have 

AA=(1-C7 m 2 )e- 2T/T2 (23) 

if 7 C m T» 1. 

It should be noted that the leading part of the out- 
put pulse can be generated without distortion only if the 
minimum value of C (t) at the initial moment of time (far 
from resonance) is larger than the lowest possible value 
of the cooperative parameter, C m /!F 2 . It follows from 
Eqs. (19) and (20) that in the case of C m > 1 we need 



T 2 > 2exp(2C m 7T). 



(24) 



This condition leads to extremely high values of the fi- 
nesse for large delay times. Fortunately, as we shall see 
below, this is no problem in a real case of storage and 
retrieval during finite time intervals. 

As noted in Introduction, the optimal strategy for 
maximizing the storage and retrieval efficiency in a quan- 
tum memory device involves the storage process that is 
the time-reverse of retrieval. Taking this into account we 
consider now the storage of the double-sided exponential 
single-photon pulse which is the time-reversed replica of 
the output pulse Eq. (17): 



(25) 



F in (t) = F ont (-t)/VJV 
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Thus, the input pulse begins at t = — oo, terminates at 
t = 0, and its maximum falls at t = —T. During the in- 
teraction between the input pulse and the atomic system 
the time-dependencies of p(t), (3(t) and, consequently, 
C(t) are assumed to be time-reversed replicas of those at 
the emission stage. In this case, from Eq. (15), assuming 
P{— oo) = 0, we obtain 

P(0) = - f F in (T)F out (-T)dr 

J — OO 

= --^L J F out (-r) 2 dr = -VAf. (26) 

The total efficiency of a quantum memory device — i.e., 
the efficiency of storage followed by retrieval — may be 
defined as 



£ = 



Jo 



(27) 



The denominator is equal to unity by definition (25), 
whereas the numerator is equal to P(Q) 2 Af so that we 
obtain £ = J\f 2 . In the limit of C m ^> 1, the efficiency be- 
comes £ w exp(— 2T/T 2 ), i.e., the losses are entirely due 
to the phase relaxation on the resonant transition during 
the delay time between the input and output pulses, IT . 

In a more realistic case, the processes of storage and 
retrieval are finite in time. Since the present scheme is 
aimed at the processing of symmetric pulses, we proceed 
by considering truncated storage and retrieval during the 
time windows [— 2T, 0] and [0, 2T], respectively, provided 
that the input pulse has the same form as Eq. (25). In 
this case 
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Whereas the norm of the infinite output pulse, Af, de- 
creases with T monotonically, the norm of the truncated 
pulse, A/trunc, which is equal to zero when T = 0, has 
a maximum at some value of T. The reason is that the 
larger T, the smaller the maximum amplitude F m , but, 
on the other hand, the smaller a truncated part of the 
pulse. It is follows from Eq. (28) that in the limit of 
C m 3> 1 , the maximum value of the norm 



(Atrunc)max ~ 1 



m(C m /2) + 1 



is achieved when the time window is 
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Thus, when the double-sided exponential pulses are ab- 
sorbed and emitted during the optimal finite time inter- 
vals, the total efficiency of storage followed by retrieval 
becomes 



£ = l-2/C m -21n(C m /2)/C m , 



(31) 



which is less than the highest possible value (C m /(1 + 
C m )) 2 ~ 1 — 2/C m . The difference is obviously due to the 
phase relaxation on the resonant transition, which in our 
model is entirely connected with spontaneous emission 
into free-space modes, during the storage and retrieval 
of information. For a given efficiency, the delay 2T max 
can be considered as a minimum (intrinsic) storage time. 
Additional storage of a photon for some time, say T s , 
simply introduces e~ Ts ^ T2 decay into the amplitude of 
the output pulse. 

In the case of the optimal time window, 2T max , the con- 
dition (24) reduces to T 2 > C m or T > aL/2n, which is 
always valid as optically thin materials are used. On the 
other hand, if T 3> aL/2ir, then only small deviations of 
p from the exact resonance — i.e., Sp/X <~ s/Cm/F <C 1 
— are necessary for implementing storage and retrieval. 

An important question is the influence of time jitter — 
i.e., the time uncertainty of receiving the input photon 
within the fixed storage window — on the efficiency. If 
we assume that the maximum of the double-sided expo- 
nential input pulse (25) falls at the time — T + S, then 
the amplitude of the excited atomic state created at the 
end of the storage process, -P(O), scales as 



P(0) oc 



r-2T 

Jo ' 



-7(C m +l)|r- 



T\ e - 7 (C m + l)\r-T- 



S Ut. (32) 



Taking the optimal storage window, T 



Tmax, we obtain 



PsMO) _ 1 + 7Cm|<5| - (1 + e 2 ^l*l) /C„ 
Pfco(0) ~~ ~ " " 



-lC m \f>\ 



1 - 2/C m 

(33) 

for C m ^> 1 and \S\ < T. According to this formula, 
the total efficiency, which scales as (P < 5^o(0)/-P<5=o(0)) 2 , 
is reduced by 10% (1%) when jSC m w 0.35 (0.1) and 
C m > 100. In other words, quantum memory efficiency of 
90% (99%) requires temporal synchronization to within 
35% (10%) of the pulsewidth. The same limitations arise 
when considering fidelity of quantum gates [37] . 

Finally, we discuss possible experimental conditions 
under which the storage and retrieval of the double-sided 
exponential pulses is possible. As a promising material 
for the quantum storage we consider rare-earth-metal- 
doped crystals (see the review in [38]) in which phase 
relaxation time at liquid helium temperature can reach 
several ms [39]. Though the inhomogeneous broaden- 
ing of resonant lines in such crystals are usually several 
GHz, a narrow absorption line on a nonabsorbing back- 
ground may be prepared using hole-burning techniques 
[40-45]. The absorption linewidth is determined by the 
laser linewidth used for the preparation and is bounded 
from below only by the homogeneous broadening. Tak- 
ing the absorption linewidth, T = 10 kHz, we obtain 
T 2 = I/71T = 32 fis. The finesse of the cavity with an 
impurity crystal is limited mainly by absorption losses 
in the host material which are typically of the order of 
10~ 2 cm -1 . So, taking a thin crystal, L = 1 mm, it 
is possible to reach T ~ 10 3 and obtain C m = 100 for 
aL ~ 10 _1 . In such a case, the duration of input and 
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output pulses is equal to 320 ns, the optimal duration of 
storage and retrieval processes 2T max w 1.25 /is, and the 
overall efficiency £ w 90%. 



IV. CONCLUSION 

It is shown that double-sided exponential single- 
photon wave packets resonant with a transition in an 
extended ensemble of atoms enclosed in a tunable cav- 
ity may be stored and recall with a near unit probability 
using cavity tuning only. The developed scheme readily 
can be combined with that for passive creation of orthog- 
onal subradiant states considered in [12], which in princi- 
ple allows one to organize multichannel quantum storage 
on a homogeneously broadened resonant transition for 
time-symmetric pulses. Besides, storage and retrieval on 
several frequencies can also be performed using periodic 
absorption structures proposed in [21], provided that the 
cavity linewidth is sufficiently large. 
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APPENDIX 

In order to obtain the solution for F(t) we integrate 
Eq. (11), which gives 

h.s(t) = /k, s (0) - i<fc, a >/j\ty(ko - k) 

x f yfZ^~T)P(T)e-^ -^ T dT. (A.l) 
Jo 

Then we substitute this result into the general formula 



F(t) 



Sc 

V 



1/2 



!>.«(*) e 



i(cJo —oj)t 



convert the sum to an integral, 
V 



E 



(2tt)3 c 3 J q 



Lo'dw / dft k , (A.3) 



and take into account that in the Born-Markov approx- 
imation, 
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(A.4) 



where the imaginary part may be neglected [46] , and 



J dft k £(d • £k, s )0(k O - k) v^CkTI) 
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lidy/CQ), (A.5) 



where /x is a geometrical factor [34] . The latter is equal to 
3\ 2 (8nS)~ 1 for a pencil-shaped cylindrical volume with 
the cross-section S, the Fresnel number F > 1, a vec- 
tor ko lying along the axis of the cylinder and dipole 
moments oriented perpendicular to the axis. In calculat- 
ing (A.5) we take into account that the diffraction an- 
gle of collective emission, which is given by the function 
(j>(ko — k) and of the order of 47r/x, is smaller than a cav . 
As a result, we obtain Eq. (14). 

To obtain the solution for P{t) we substitute the solu- 
tion Eq. (A.l) into Eq. (12), employ Eq. (A.3) and take 
into account that 



J dn^id ■ £k,,) V(k - k)£ k (t) 



8?r 



(H+ l/N)d 2 C{t). (A.6) 



This leads to the following equation 

^ = - 7 (l + C(t))P(i)-^)F in (t), (A.7) 



(A. 2) which has a general solution of the form (15). 
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